Quarkonia Suppression in QGP Under Colour Screening Scenario at SPS, 

RHIC and LHC 



3 



P. K. Srivastava a , M. Mishra h and C. P. Singh a 

° Department of Physics, Banaras Hindu University, Varanasi - 221005, India 
b Birla Institute of Technology and Science, Pilani - 333031, India 

Abstract 

We present a modified colour screening model for J/ij) suppression in the Quark-Gluon Plasma (QGP) 
using quasi-particle model (QPM) as equation of state (EOS). Other theoretical ingredients incorporated 
in the model are feed-down from higher resonances namely, \ c , and ip , dilated formation time for 
quarkonia and viscous effects of the QGP medium. Assuming further that the QGP is expanding with 
Bjorken's hydrodynamical expansion, the present model is used to analyze the centrality dependence of 
the J/ip suppression in mid-rapidity region and compare it with the data obtained from SPS, RHIC and 
LHC experiments. We find that the centrality dependence of the data for the survival probability at all 
energies is well reproduced by our model. We further compare our model predictions with the results 
obtained from the bag model EOS for QGP which has usually been used earlier in all such calculations 
and our results indicate that the present EOS yields a better result in comparison to the bag model. 
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1 Introduction 



Ultra-relativistic collisions of heavy nuclei are believed to produce nuclear matter at extreme pressure and 
energy density. Quantum Chromodynamics (QCD) predicts that a deconhned state of partonic matter which 
is referred as the quark-gluon plasma (QGP), will be produced in such circumstances. One of the most 
striking signatures is the suppression of heavy quarkonia states such as the charmonia {J/ip, ip , Xc etc) in 
the deconhned state. Charmonia suppressions in heavy-ion collisions have been experimentally studied, first 
at the CERN Super Proton Synchrotron (SPS) by NA50 [2] and NA60 experiment [3] at ^sWn = 17.3-19.3 
GeV and then at Relativistic Heavy Ion Collider (RHIC) by the PHENIX experiment at v / J/vw = 200, 39 and 
62.4 GeV g], and now at Large Hadron Collider (LHC) by the CMS and ALICE experiments at ^/snn = 2.76 
TeV [SJIS]- However, experimental results involve more puzzling features which defy explanations on the 
basis of colour-screening in the QGP alone. For example, we notice significantly less suppression at the 
mid-rapidity than at forward rapidity (1.2 < \y\ < 2.2) while the medium in the central part of the collision 
is the most dense and intutively one sould expect an opposite behaviour. Another interesting experimental 
observation is similar J/ip suppression at SPS and RHIC energies for same number of participants 0]. 
Many theoretical efforts have since been made which explore the effects of cold nulear matter (CNM) as 
well as possible coalescece of originally uncorrelerated c , c quark pairs. However, the complete theoretical 
interpretation for J / ip suppression is still lacking. 

Many theoretical papers have analyzed the heavy-ion experimental data regarding nuclear modification factor 
for heavy quarkonia. For example, Pb + Pb data at SPS explained by variety of models [7J [5J HI OH EJ 
[T2J [T31 H3] with or without deconfinement phase transition scenario. Yunpeng et al., [TS] analyzed the J/ip 
suppression data using a transport model calculation involving suppression as well as the regeneration of 
J I ip due to coalescence of c and c pair produced abundantly in the initial stage of collisions at RHIC and 
LHC energies. Zhen Qu et al., |16) used a similar approach to explain the forward rapidity data for J/ip at 
RHIC energy. Sharma and Vitev [17] calculated the yields of quarkonia at RHIC and LHC as a function 
of transverse momentum based on heavy quark effective theory in which both colour singlet and colour 
octet contribituions along with feed down effects from excited states were suitably incorporated. Their 
model provides good description of the LHC data for central as well as peripheral collisions. But none 
of the available models could provide a satisfactory and consistent explanation for the suppression at all 
energies simultaneously. We have recently modified the colour-screening model [TB] of Chu and Matsui 
by parametrizing pressure instead of energy density since pressure density becomes almost zero at the 
deconfinment phase transition point. We have further shown that the variation of survival probability of 



J/tp with respect to participant number obtained from this model compares well with the experimental 
data (CNM normalized) at RHIC [IS]. In the present paper, we reinvestigate the colour-screening model by 
employing a more realistic thermodynamically-consistent quasi-particle prescription for the equation of state 
(EOS) of QGP medium. This highlights a major difference between our present approach and the model 
of Chu, Matsui [Jj|] or Mishra et al [TB] because bag model is still considered as a crude EOS for QGP. We 
have demonstrated earlier that various thermodynamical as well as transport quantities obtained from the 
present quasiparticle model (QPM) compare well with the most recent lattice results [30] [H|. In addition, 
we incorporate here viscous hydrodynamics in our formulation with a constraint that the ratio r\js varies 
very slowly |21) with the temperature in the QGP medium. This assumption helps us in taking r\ to be 
almost independent of the temperature. Feed down from higher resonances have also been incorporated in 
the model. 

We use our formulation to explain the recent LHC experimental data on J/tp suppression along with the 
RHIC and SPS data. We surprisingly find that J/ip suppression at mid-rapidity is explained consistentely 
by our model at SPS, RHIC and LHC energies in QGP picture alone. Further we compare results obtained 
from quasiparticle model (QPM) and the bag model EOS on J/ip suppression and find that QPM EOS 
describes the data better in comparison to the bag model EOS for QGP and thus supports the utility of the 
present model in explaining J/ip suppression. 



2 Formulation 

2.1 Cooling law 

The main theme of the present formulation is essentially borrowed from our previous model [18] where we 
have used the colour screening idea of Chu and Matsui [TH]- We have assumed that QGP medium formed 
during the collision, expands and cools according to the Bjorken's boost invariant longitudinal hydrodynamics 
in mid-rapidity region. Employing the conservation of energy-momentum tensor, the rate for the decrease 
of energy density e 22 is given by 

de e+p irj 

d V t 3r 2 ' lj 

where r\ is the shear viscosity of the QGP medium, p is the pressure density and r is the proper time. Using 
Eq.(l) and the thermodynamical identity e = — p, the cooling laws for energy density and pressure can 



be separately derived as (see Appendix) 
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where c\ , ci and C3 are constants which can be determined by imposing the initial boundary conditions on 
energy density and pressure, q = c 2 s + 1 with c s as the speed of sound in the medium. The energy density 
and pressure are computed by employing quasiparticle model (QPM) EOS for QGP. We take e = eo at 
t = To (initial thermalization time) and e = at r = r ; where r is the proper time. Consequently, the 
constants C\ and ci are related as : 

, where r = tqA^ 1 ^ 1 , A = T /T and i? is the Reynold's number for QGP. Further : 
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Using the other initial condition p = po at r = ro gives the value of C3 as follows : 



/ , s £ 2-1 q \ (ci-i) 



2.2 Pressure Profile 

Assuming that the pressure almost vanishes at the transition temperature T = T c , we take a pressure profile 
function in the transverse plane with a transverse distance r as }18j : 



p(U,r)=p(U,0)h(r); h(r) = ( 1 - ^ ) 9(R T - r), (7) 



where the coefficient p(ti, 0) is yet to be determined, Rt denotes the radius of the cylinderical plasma and it 
is related to the transverse overlap area At as determined by Glauber model Rt = y^F' ^ e pressure is 
thus assumed to be maximum at the central axis but it vanishes at the edge Rt where hadronization occurs. 
The exponent /? depends on the energy deposition mechanism and here we have taken ft — 1.0 |18j ; 9 is the 



unit step function. The factor p(U, 0) is related to the average initial pressure < p >j via 



p(t i ,0) = (l + l3)<p> i . (8) 

The above average pressure is determined by the centrality dependent initial average energy density < e >j 
which is given by Bjorken's formula [231 121] : 

1 dE T 

< e >i= j- , 9 

Am dy 

where dET / dy is the transverse energy deposited per unit rapidity. We use the experimental value of dET / drj 
where r\ is the pseudorapidity and then multiply it by a corresponding jacobian factor |24[ 125) in order to 
obtain dET/dy for a given number of participants {N part ) at a particular center-of-mass energy (y^s^T/v)- 
Further at the initial proper time fff^- = ^ L = c 2 s in QPM EOS [2T] and thus < p >;= c 2 s < e >,-. 

2.3 Constant Pressure Contour and Radius of Screening Region 

Since the cooling law for pressure cannot be solved for r analytically and therefore we use a trick to determine 
the radius of screening region. Writing the cooling law of pressure as follows : 

/ x , BCD 
p(T,r)=A+ — + - + —, (10) 

where A, B, C and D are constants related to c%, C2 and C3 as : A = — c±, B = C2C 2 ., C = ^2!^ and 
D = C3. Writing the above equation at the initial time r = Ti and at screening time r = r s we get : 

BCD 

p(n,r)= A + — + - + — ^=p(n,0)h(r), (11) 

and 

, BCD 

p(r s ,r) =A+ — + - + —= p QGP . (12) 

Ts T s Ts b 

Here p QGP is the QGP pressure as determined by EOS in QPM [20l [21] . Solving Eqs. (11) and (12) 
numerically and equating the screening time t s to the dilated formation time of quarkonia tp {—^Tp where 
7 = Et/M^ is the Lorentz factor associated with the transverse motion of the c — c pair and = 3.1 GeV 
and tf is the proper time required for c — c for formaition of J/^) , we can find the radius of the screening 
region r s . The screening region is defined as a region where temperature is more than the dissociation 



temperature so that the quarkonia formation is unlikely inside that region |18) . Hence the pair will escape 
the screening region and form quarkonia if \r^ + vtp\ > f s where is the position vector at which the 
charm-quark pair created. 

The above kinematic condition takes a simplified form by assuming that J/tj) is moving with transverse 
momentum px- Thus the above escape condition can be expressed as a trigonometric condition [18] : 

cos(j)>Y; Y= 1 , (13) 

where cj> is the angle between the transverse momentum (j>t) and the position vector f\p and ru = \rJ\ with 
m = M$. 



2.4 Survival Probability 



Assuming the radial probability distribution for the production of cc pair in hard collisions at transverse 
distance r is 

(1-^r) 0(R T ~r). (14) 



Here we take a = 1/2 in our calculation as used in Ref. |19j . Then, in the colour screening scenario, the 
survival probability for the quarkonia can easily be obtained as [18j [19] : 
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where the maximum positive angle 4> max allowed by Eq. (13) becomes 



(15) 
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Generally experimentalist measure the quantity namely, px integrated nuclear modification factor and there- 
fore, the theoretical pr integrated survival probability is of paramount importance to compare with the 
experimental results. It is defined as follows : 



s'(n , ^::: s '^T^ N part)dPT 

O {Apart) - j-PT^ax ^ • U b j 

It has been found that only about 60% of the observed J/ ip come from hard collisions whereas 30% is from 
the decay of \c and 10% from the ip [IB]- Therefore, the net survival probability of J/ip in the presence of 
QGP medium is : 



S ifi> = 0.605^ + 0.303^ + O.IOS^ (17) 

In our calculations we use T c = 0.17 GeV in accord with the recent lattice QCD results, f3 — 1.0, and 
a = 0.5 [IB]. We take the masses as 3.1, 3.5 and 3.7 GeV of J/ip, Xc and ip , respectively [2E]- Further we 
choose formation time tf as 0.89, 2.0 and 1.5 fm for J/ip, Xc and if) , respectively [55J[TB]. We use these 
values of parameters at all y/SNN m our calculations for survival probability. 
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Figure 1: Variation of survival probability of J/tp at RHIC energy {Sf^ IC ) with respect to N part . Experi- 
mental data points are taken from Ref. [4]. 



3 Results and Discussions 



Fig. 1 shows the variation of survival probability of J/ip at RHIC energy (i.e. -Jsnn = 200GeV^) with 
respect to N part . Star symbols represent the results of Ref. [18] where a bag model EOS for QGP along 
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Figure 2: Variation of survival probability of J/ip at SPS energy {SjJ^ s ) with respect to N part . Model I and 

Model II are our model with two different sets of dissociation temperatures (Tj) for J/tp, Xc and ip (details 
in text). Experimental data points are taken from Ref. [2]. 
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Figure 3: Variation of survival probability of J/tp at RHIC energy (S^, IC ) with respect to N part . Experi- 
mental data points are taken from Ref. [4]. 



with ideal hydrodynamics was used and S^R 1 has been calculated at fixed pr = 2.0 GeV with dissociation 
temperatures Td/T c — 2.1, 1.16, 1.12 for J/tp, Xc and ip , respectively [18]. Open circles are our model 
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Figure 4: Variation of survival probability of J/ip at LHC energy (Sf/^ C ) with respect to N part . Experimental 
data points are R e ^A f° r prompt J/ip without any CNM normalization as taken from Ref. [5] . 
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Figure 5: Variation of survival probability of J/ip at RHIC energy {S^ IC ) with respect to N part calculated 
within our model in viscous as well as in ideal hydrodynamics framework. 



results where we have used QPM EOS along with viscous hydrodynamics and parameters have the same 
value [18 . The comparison shows that the results from our model provide a better description of the 
experimental data [IJ (shown by solid points in Fig. 1) in comparison to the results obtained earlier |18j . 



Here experimental data have been normalized with CNM effect. We have also shown the pr integrated 
variations of S^^ IC with N part obtained from our model, represented by solid inverse triangles, using the 
same values of dissociation temperature and observe that the agreement with the experimental data has 
improved. 

Fig. 2 presents the variation of SjEf with N part obtained from our model for two different sets of dissociation 
temperatures (labelled as Model I and Model II ) and their comparison to the CNM normalized experimental 
data [HUB]. Moreover, one should keep in mind that from here onwards, we use the same px range in all of 
our model calculations as used in related experimental data. In model I (shown by filled square), we have 
taken the dissociation temperatures Td/T c — 2.1, 1.16 and 1.12 for J/ip, Xc and -0 [IB] whereas model II 
(shown by solid triangle) gives the calculation for T c i/T c = 2.0, 1.5 and 1.4 for J/tp, Xc and V' > respectively. 
In some recent lattice calculations [26, 27], after constant mode correction, there is no hint for suppression 
of Xc and ip instantly after T c as pointed out by earlier studies 28, 29, [301 13H E2 ■ Due to above uncertainty 
regarding the dissociation temperatures of different quarkonia states, we wish to see the effect of Td in 
explaining the data of RHIC and SPS experiments, simultaneously. We observe that the changes in Td 
values produce more significant results in the peripheral collisions rather than in central collisions. 

Fig. 3 demonstrates the variation of S^K 1 with respect to N part as obtained from model I and model II. 
Here we observe that except for extreme peripheral collisions, model II compares well with the experimental 
data (normalized with CNM effect) in comparison to model I for semi-central case and results from both 
models converge for central collisions. This again shows that a change in Td values is more effective in 
peripheral collision rather than in the central collision. However, at RHIC energy the difference in the 
results of two models is insignificantly less in comparison to their results at SPS energy for all collision 
centrality. It may be due to the reason that when we increase the dissociation temperatures for Xc and ip' , 
the energy deposited in the collision at SPS energy is still not sufficient to melt all the XcS and i/j's even in 
the central collisions and thus it provides a significant contribution to the survived J/ips. However, at RHIC, 
except for most peripheral collisions, the energy deposited in the collision zone becomes always sufficient 
to dissociate a significant fraction of XcS and ip's even if their dissociation temperature becomes higher as 
1.5 T c and 1.4 T c from 1.16 T c and 1.12 T c , respectively. Thus the contribution from Xc and ip' is less in the 
survival of J /ip. Moreover, for most central collisions the contributions from Xc and ip' becomes negligible 
and the prediction of both the models overlap on each other. 

Fig. 4 shows the variation of S^^ c with respect to N part obtained in our model I and II. Here we find 
that there is almost no difference in the survival probability calculated within model I and model II for all 
collision centrality. This means that the contribution to Sj/^ from % c and ip' is almost negligible and the only 



contribution to survival probability comes from J / ijj alone. We also present the recent experimental data of 
nuclear modification factor (Raa) without any cold nuclear matter (CNM) effect correction for prompt J/ip 
in mid-rapidity (|y| < 2.4) obtained by CMS collaboration 5.. The qualitative agreement of our prediction 
with the experimental data is reasonably good and quantitative agreement is expected only after normalizing 
the experimental data (nuclear modification factor) by the corresponding contribution from CNM effect (i.e., 
R AA M ) at LHC energy. 

In Fig. 5, we present our model results for S^/ c with N part and also show a comparison with results 
obtained in our model when we put shear viscosity (rj) equal to zero i.e., considering ideal hydrodynamic 
framework. We find that if we choose the value of rj/s = 1/A-k (suggested as the lower limit for any fluid [33 ) 
then the difference between ideal and viscous case vanishes. However, for larger values of ij/s, the effect is 
expected tobe sizable. We have also checked the difference between ideal and viscous case at SPS energy 
and find that the effect is further reduced in comparison to the results at the RHIC energy. 

In summary, we have presented a modified colour-screening model for J/ip suppression in QGP medium 
where quasi-particle model equation of state for QGP, feed down from higher resonance states (namely, Xc, 
and ip ), dilated formation time for quarkonia and viscous effect in the QGP medium have been properly 
taken into account. Furthermore, we assume that the QGP is expanding with Bjorken's boost invariant 
longitudinal expansion. The above model has been employed to analyze the centrality dependence of the 
J/ip suppression data from SPS, RHIC and LHC experiments. We find that the centrality dependent survival 
probability is adequately reproduced by the current model in a consistent manner. We have also compared 
our present model results with the presious model results (TBJ where bag model EOS of QGP has been used. 
We notice that present model describes the experimental data better as compared to the old model [IS]. 
Thus we conclude that our model provides a better explanation for the J/ip suppression at SPS, RHIC and 
LHC experiments without any free parameter. Further extension of this work is underway to explain the 
complete rapidity dependence of the data in order to draw a firm conclusion regarding the overall level of 
suppression arising from QGP picture alone. 
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Appendix 



(l+l)-Dimension Cooling Law for Temperature 

To detrmine the explicit proper time (r) dependence of temperature (T), we first start with the entropy 
equation for (l+l)-dimension expansion in (3+l)-dimension space given as follows [PRD 41, 2903 (1990)]: 



where s is the entropy density of the system and R is the Reynold's number defined as follows for this system 
having shear and bulk viscosity as r\ and (, respectively: 



R 



-l 



4/3 7? + C 

T ST 

4 V , 

3T77 (? 



WC = o). 



(19) 



After differenting Eq. (18) with respect to r we get : 




(20) 



We can rewrite the above Eq. (20) as follows : 




(21) 



After integrating both side of the Eq.(21), we get : 




(22) 



or 




(23) 



Using the boundary condition : T = Tq at r = to in Eq. (23), we get the value of C as : 



Finally, using the above value of constant C in Eq. (23), we get the 1+1 cooling law for temperature as 
below : 

/ T \ 3 / r\ (i-V-R) 

(25) 
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(l+l)-Dimension Cooling Law for Energy Density 

The rate equation for energy dansity in Bjorkcn longitudinal expansion scenario for viscous case is as follows 

de £+p Ar > (oa\ 

d~T=-— + ^ (26) 

and the thermodynamical relation between e and p in quasiparticle model (QPM) can be given as below [19]: 
Differentiating (27) w.r.t. r, we get the following relation : 

d 1 = drdP d^p dp 

dr dr dr dr 2 ' dT { ' 

After equating Eq. (28) with Eq. (26) and putting e + p = T^Jp, we got the following equation : 

_Tdp 4 !L= cPp dp ) dT_dp 

t dT 3r 2 1 dT 2 dT' dr dr ( ' 



Writing Eq. (26) as : 



de T dp Ar, 



Differentiating the above Eq. (30) two times and eliminated (T^^ + f|0§7 with the help of Eq. (29), we 
get 



d 2 e de Art dp „ ,„„. 



We can rewrite Eq. (31) as : 



Solution of Eq. (32) is given by : 

3c 2 T 



e = c 1+ c 2 r-(^) + ^i (33) 



(l+l)-Dimension Cooling Law for Pressure 



Again start from energy density rate equation : 



It can be further written as : 

de dp _ _ e _ p Aij 

dpdr~ r r + 3r2- [6b) 
Since in QPM EOS, we know that dp/de — c^, thus the above Eq. (34) can be rewritten as : 

< 36 » 

Now, using the value of e from Eq. (32), we get the final differential equation for the cooling law of pressure 
is as follows : 

To solve this above partial differential equation, we need to determine its integrating factor (I.F.) which 
comes out to be r c =. Further we can write the solution of Eq. (36), by using this I.F. as follows : 



c l Cl + C 1 C2T -^) -p^d + cl)) T ot dT (38 ) 
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P = -ci + c 2 > +-LI -|— +C3T- C ». (39) 
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P=-ci + c 2 4 + ?f^-r) +c 3 r- c », (i/ ? = ^ + l) (40) 
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Evaluation of constants cl, c2 and c3 



Using the first boundary condition e = e at t = t in Eq. (33), we get the first constraint relation as follows 

An 1 

eo = Ci+c 2 T 9 + -4- (41) 
3c| r 



Similarly, using second boundary condition e = at r = r (corresponds to temperature T = 160MeV) in 
Eq. (33), we get the second constraint relation as follows : 



= ci + c 2 
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(42) 



Subtracting Eq. (42) from Eq. (41), the expression to calculate c 2 is given below : 
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where r is calculated from Eq. (25). Consequently c\ can be calculated from the following relation : 



ci = -c 2 r' 
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Finally, using initial condition p = p at r = t in Eq. (40) , we can get the following expression to calculate 
c 3 as : 

/ | \ c i 2 -i 4? 7 ( Q \ (c 2 s -i) 
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